Pairing Interaction in N=Z Nuclei 
with Half-filled High-j Shell 



A.Juodagalvis* 

Mathematical Physics Division, Lund Institute of Technology, 
S-22100 Lund, Sweden 



Abstract 

The role of L=0 pairing interactions (both T=0 and T=l) in three selected 
bands with the total isospin T tot =0 in 22 Na, 48 Cr and 90 Rh nuclei is discussed in 
the spherical shell model. These bands were selected requiring termination in a 
most unfavourable way. The investigated effects coming from pairing interactions 
include the change of nuclear wavefunction, quadrupole properties and state 
energy. To follow the gradual change of nuclear properties, the investigation is 
performed in a perturbative approach by defining the Hamiltonian with a variable 
content of the pairing interaction. It is shown that the pairing interaction does 
not affect the quadrupole properties, if the energy distance to other close lying 
states is large. The calculated pairing energy is shown to have a similar spin 
dependence in these three bands. In addition, the influence of the model space 
as well as of the residual interaction on the derived pairing energy is discussed. 

Key words: Shell model; Pairing energy; Schematic pairing interaction; Band termi- 
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1 Introduction 



The pairing force plays an important role for nuclear properties like binding energy, de- 
formation, moment of inertia. It is a long standing question what is the quantitative 
role of the isoscalar (T=0) and isovector (T=l) pairing interactions. Their competition 
may be observed in the double-odd N=Z nuclei, where the lowest states with the total 
isospin T tot =0 and 1 have very similar binding energy, and the transition from the T tot =0 
to T to t=l ground state, as the mass number increases, is seen. This phenomenon has 
recently regained attention, see e.g. [|l| and references herein. 

To investigate the role of the pairing interactions, many attempts were made to incor- 
porate those correlations in the single-particle models like Hartree-Fock (HF) approxima- 
tion (see e.g. @, ^|, ||, U). In the HF approach, the pairing correlations should be added. 
A different situation occurs in the shell model which takes into account the residual inter- 
action completely. Here the pairing force is a part of the full interaction, and in principle 
it cannot be separated from the latter. It is known that the shell model Hamiltonian, 
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consisting of the single particle term and the T—l pairing and quadrupole interactions ||], 
gives a reasonable description of nuclear properties. However, a question whether this is 
a sufficient minimal shell model Hamiltonian appeared again, when the theoretical back- 
ground for strengths of pairing and quadrupole forces was reconsidered and it was shown 
that a schematic T=0 pairing is also an important part of the realistic interaction [[/J. 
The interest was also stimulated by the availability of advanced shell model calculations 
(see e.g. [|], ||) as well as by new experimental data available (e.g. ||, jlQl). 

This paper discusses the role of a schematic L=0 pairing. The two-particle states with 
spin S and isospin T equal to ST=01 or 10 are the ones favoured by Wigner's SU(4) 
classification. Those states are also among the few lowest eigenvalues of the multipole 
part of a realistic Hamiltonian J7j. The real nuclear interaction breaks SU(4) symmetry 
(e.g. by spin-dependence). However, it is not completely destroyed, and its presence shows 
up most clearly in the lighter nuclei through selection rules for /3-decay. 

The influence of the isoscalar and isovector L=0 pairing components of the effective 
interaction has been discussed for the nuclei in the middle of p/-shell |L[ || lT| . The 



present paper refers to three N=Z nuclei having half-filled high-j shell. These three 
selected nuclei ( 22 Na, 48 Cr and 90 Rh) are of particular interest. First of all, since protons 
and neutrons are occupying same shells, the neutron-proton interaction plays an important 
role. Another interesting feature is a rotational behaviour of the bands that terminate at 
spin I m =(p+l) 2 h (p is the principle quantum number of the harmonic oscillator shell). 
Those bands show backbending in the spin region (J m — 4,J m — 2). They terminate in a 
spherical state, since all components of the high-j shell orbits having positive angular 
momentum projection are occupied. In addition, they terminate in a "most unfavoured" 



way |L2|, since last two transitions before the termination involve breaking and alignment 
of two T—l pairs (one neutron pair and one proton pair). Besides these three nuclei, such 
a band termination may occur in other nuclei (for example, an excited positive parity band 
in 47 Cr, having one hole in the sd-shell), but those bands would lie at higher excitation 
energy. 

The most unfavoured band termination was the criterion to select particular bands in 
these three nuclei. These bands have the total isospin T tot =0. The role of pairing interac- 
tions was investigated using the shell model. For two nuclei, 22 Na and 48 Cr, experimental 
data exists. But no data is available for 90 Rh, thus we have to rely on the ability of the 
chosen interactions to describe properties of the neighbouring nuclei. The aim of this 
paper is to show that the selected bands have some features in common, and that the 
L=0 pairing interactions change the properties of the bands in a similar way. The pairing 
influence to the nuclear wavefunction and quadrupole properties is also discussed. 

The paper is organized in the following way. Computational details and comparison 
with experimental data are both discussed in section |2|. Similarities in the selected bands 
as well as the predicted properties of the selected band in 90 Rh are also discussed in a 
greater detail there. At the end of the section, evaluation procedure of the role of pairing 
interaction is explained. The results are discussed in section |3|. It is shown that the 
L=0 pairing interaction affects the behaviour of the selected bands in a similar way. And 
finally, conclusions are drawn in section 

2 Description of the approach 

This section discusses three selected bands in a greater detail (subsection |2.1| ). Then 
computational details and comparison with experimental data are given (subsection |2.2|). 
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Finally, an approach to extract the role of pairing interaction is explained (subsection 
|2~3| ). The obtained results are discussed later in section |3|. 

2.1 Bands having most unfavoured termination 

Among the N=Z nuclei, there is a special class of nuclei having exactly half-filled high-j 
shell. A peculiarity in these nuclei is that it costs more energy to construct the state, where 
all the valence particles are aligned in this high-j shell giving maximum spin I m =(p+1) 2 
(where p is the principle quantum number), than a rotational trend would suggest [fL^] . 
This effect comes mostly from JT =01 pairing interaction. It is interesting to see that 
there are similarities among the bands also in the spin dependence of the content of 
pairing energy. 

In this paper, three N=Z nuclei are discussed, namely, 22 Na, having 6 valence particles 
(3 protons and 3 neutrons) in the d 5 / 2 , 48 Cr, having 8 particles in / 7 / 2 , and 90 Rh, having 
10 particles in gg/2 shell, when the single-particle level ordering from a harmonic oscillator 
potential with £ 2 and t ■ s terms is assumed. In the single- j shell model with the JT=01 
pairing interaction only, bands leading to the maximum spin state behave in a similar way 
(see fig.|]a). Let J m is the spin of a completely aligned state (which should be spherical 
since all orbits with positive j projections are occupied): I m =9h for 22 Na, 16h for 48 Cr, 
and 25h for 90 Rh. Let us discuss only the lowest states at spins I m —2k where k is an 
integer number. These states have the total isospin T tot =0. The state with spin J m — 2 
would be a linear combination of states where one pair of identical particles is coupled 
to angular momentum 0, and the spin is obtained from the remaining unpaired particles. 
Similarly, the J m — 4 state would have two pairs of particles (one neutron pair and one 
proton pair) coupled to angular momentum 0. For 22 Na, that is the maximum spin state 
produced by two odd particles (5 + , bandhead K n =3 + ). Since there is no change in the 
pairing energy when going from J m — 4 to I m — 6 state, the I m — 4 state may be somewhat 
irregular in the real energy spectrum. In the discussed cases there is a backbending in 
48 Cr and 90 Rh, but in general the change of the moment of inertia may be less pronounced. 
The state with spin I m — 10 (having 10 units of JT=01 pairing energy in a single- j model, 
see figfja) would be a linear combination of states where 6 particles contribute angular 
momentum zero: two pairs of one kind and one pair of the other kind of particles are all 
contributing angular momentum 0. Since at least 6 valence particles are not contributing 
to spin, the irregularity in the energy spectrum may be not seen unless this j-shell is 
dominating. The state J m — 16 would have 4 angular-momentum-zero pairs. This state is 
the ground state in 48 Cr, while in 90 Rh it is the maximum spin obtained from odd proton 
and neutron (9 + , bandhead K n =5 + ). 

There are more bands terminating in a most unfavourable way. One can think of lighter 
nuclei with half-filled high-j shell of harmonic oscillator major shells having principle 
quantum number p=l or ( 8 Be and 2 H, respectively), but there are very little states 
before the maximum spin state, thus no bands are formed. A more realistic examples 
are excited bands in nuclei close the three earlier mentioned nuclei. First example could 
be one-particle-one-hole (lp-lh) bands of unnatural parity in nuclei having mass number 
(A m —1) (here A m is the mass number of any of the three middle high-j shell nuclei). For 
example, an excited band in 21 Ne (or 21 Na), where one proton (neutron) is excited from 
the p- to s(i-shell. This band would terminate at spin I=9+jh, where jh equals 1/2 or 3/2 
depending on the position of the hole in p-shell. Other examples are 47 Cr (or 47 V) and 
89 Ru (the band in 89 Rh is less likely to be observed, at least in the nearest future). One 
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Figure 1: Spin dependence of the JT pairing energy in the lowest T to t=0 states in the 
model space consisting of a high-j shell only. The JT=01 case is shown on the left, while 
JT=10 case is shown on the right. Formula for the JT=01 pairing energy, parameterized 
in terms of the total seniority and reduced isospin, is taken from fl3 |. 



can also think of bands based on 2p-2h excitation in (A m —2) mass nuclei, and so on. 

The most unfavoured band termination may also be encountered in cases when the last 
transition before the band termination involves alignment of a pair (coupled to angular- 
momentum-zero) that brings only 2h, and the relative difference of neutron and proton 
seniorities (i.e. the reduced isospin) decreases. In the cases discussed in this paper, the 
reduced isospin is zero at spin J m — 4, then it increases at state J m — 2 and decreases again 
at the aligned state with the total angular momentum J m . 

The discussion of pairing features in a single- j shell is based on the fact that the high-j 
shell dominates in the cases discussed here, i.e. either it has the largest occupancy and/or 
its change has greater impact to the band properties. For example, shell occupancies in 
some states of 48 Cr are shown in fig.[7| (occupancies in all yrast states of 48 Cr may be found 
in ref.||14|). The shell occupancies along the selected band in 90 Rh are shown in ng.fl 



, see e.g. 
are less 



The role of JT=01 pairing interaction is well acknowledged. The seniority quantum 
number led to an additional classification of many-particle configurations. While the role 
of JT=10 pairing interaction is less pronounced: models without this interaction suc- 
cessfully describe nuclear properties, and there is no additional (and useful) classification 
in terms of JT=10 coupled pairs. Although in the realistic interactions JT=10 pairing 
is stronger than JT=01 0: the ratio of the strength parameters Gio/Goi~1.5 
Table [I]; the features of this pairing interaction in a single- j shell (see fig.0fc) 
prominent in the calculated or observed spectrum. For example, the behaviour of JT=10 
energy in a single- j model predicts an increased energy loss between states having spins 
I m — 6 and I m — 4. In the observed spectra (fig.||) as well as in the full calculations (figs.|],|3]), 
there is a backbending between those states, i.e. the behaviour is opposite to the predicted 
from the JT= 10 pairing energy spin- dependence in a single high-j shell. It may be noted, 
however, that the unequal energy spacing between states with spins J m — 4, I m — 2 and I m 
is predicted by JT=10 as well as by JT=01 pairing interaction. 

One reason for the small relative importance of JT= 10 pairing may be the fact, that 
the interaction matrix elements of JT =10 pairing are smaller than those of the JT=01: 
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for example, in the single- j case, the ratio of the matrix elements is G w /Gqi (j + l)/(3j) 
which is less than 1 for j>\- The pairing energy would be larger if the high-j shell 
spin-orbit partner would be included in the model space, as discussed below in subsection 



3.2.1] . However, that does not change much the spin dependence of the pairing energy (see 
figjfj]). Thus the main reason, why the JT= 10 pairing interaction was not found being 
important, may be its smooth spin dependence (see fig.[3]&). Small irregularities suggested 
by fig.|l|& are likely to be found only in the cases when JT =01 pairing energy does not 
change. Another possibility could be the cases when JT=01 pairing is blocked, e.g. in 
the doubly-odd nuclei |pT| . In subsection |3.2.2j it is discussed that the JT= 10 pairing 



interaction may be responsible for the favouring of 1=7 state in 90 Rh calculated using 
SSSV interaction. 

The discussion above was referring to the pairing energy spin dependence along one 



band, and is continued in subsection |3.2.3| . One may find a more general discussion 



about the relative importance of JT=01 and JT=10 pairing energies in the shell model 
framework in ref. 

2.2 Computational details and comparison with experiment 

In present work, shell model calculations were performed using the computer code AN- 



TOINE ||15|| . The residual interactions, used in the calculation, were: USD |T6| interaction 
in the full s<i-shell for 22 Na, KB3 [[17J in the full p/-shell for 48 Cr, Gross and Frenkel in- 
teraction [[Tj| (referred to as GF) as well as SLGTO |19[ in pi/25'9/2, and Sinatkas et al. 



[pOfl (referred to as SSSV) in (79/2P1/2P3/2 for 90 Rh. Electromagnetic quadrupole properties 
were calculated using the effective charges q n =1.5e for protons and q u =0.5e for neutrons. 

Fig. ^presents a comparison of the calculated and experimental energies of the selected 
T to t=0 states in 22 Na and 48 Cr: These are odd-spin members of the ground state band 
having K 7T =3 + in 22 Na (experimental data taken from jnj]), and the ground state band 
having K 7T =0 + in 48 Cr (experimental data taken from ||22|| ). As can be seen, the agreement 
between experimental data and calculated values is excellent. There is no experimental 
data available on 90 Rh. In fig.§ the calculated odd-spin members of the K n =5 + band 
in 90 Rh are shown. Those two figures (1 and 2) show, that after a rotational reference 
E TOt is subtracted, the bands tend to bend up close to the band termination (so-called 
"unfavoured band termination" f23|). In addition, the state I m — 4 deviates from a regular 
behaviour of the preceding states (not seen in 22 Na, since there is only one such state). 
The irregularity at spin J m — 10, discussed in subsection ^4] as being possible to occur, 
is not seen in 48 Cr and 90 Rh at 1=6 and 15, respectively. The exceptions are GF and 
SLGTO calculations of 90 Rh, where a very restrictive model space allows rather poor shell 
mixing. The calculated probability of a configuration, where all active particles are in the 
high-j shell producing spin J m , is large: >90% for 48 Cr as well as 90 Rh (SLGTO and SSSV 
calculation), 54% for 22 Na, and 62% for 90 Rh (GF). 

Since there is no experimental data on 90 Rh available, only calculated states are shown 
in fig.|3]. Calculations were performed in the configuration spaces (pi/2<?9/2) 14 with GF and 
SLGTO (cf.|fr9[) interactions, and (5 , 9/2Pi/2P3/2) -10 using SSSV interaction. These three 
residual interactions predict that the ground state has spin 1=0 and the total isospin 
T tot =l, in agreement with the mass-dependence of the both spin and isospin in the ground 
state of doubly-odd nuclei PS| . The excitation energy of the lowest T tot =0 states is 
predicted to be ~1.3 MeV. The SSSV interaction predicts almost spherical ground state 
of 90 Rh, with the deformation (3^0.13 (value estimated from B(E2;2i 0f,T tot = 1), 
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Figure 2: Energies of selected bands in 22 Na and 48 Cr. For the calculated states, the resid- 
ual interaction is given in parentheses. Calculated spectra of 22 Na and 48 Cr were reported 
in refs. fl6|] and ||24fl , respectively. Plot a shows excitation energies, and plot b shows the 
same energies but with a rotational reference, -E'rot=32.3274~ 5 ' /3 /(/+l), subtracted. In the 
latter plot, 4.5 MeV was added to all energies in 22 Na to make the figure more compact, 
that is also indicated by an arrow. 

assuming axially symmetric rotor model is valid). The interactions, however, predict 
slightly different behaviour of the band of interest. For example, the SSSV interaction 
favours I 7T =7 + state as compared to the 5 + state. In addition, energies calculated in a 
more restricted model space (with the GF and SLGTO interactions) are more irregular 
than those obtained using the SSSV interaction, see fig.|3]&. Like in 48 Cr case at spin 
I m — 10=6 + , there is no irregularity at spin I m — 10=15 + in 90 Rh according to the SSSV 
calculation. That may be explained by the shell mixing (see occupation numbers in fig.§). 
The band terminates at spin 7 m =25 + in an almost spherical shape, with a half-filled gg/ 2 
shell. An exception is the GF calculation where the terminating state contains 62% 
probability of configuration {p\/ 2 9gj 2 ) 2 an d 37% of (p\, 2 gg/ 2 ¥ and the nucleus has an 
oblate shape. The band termination is also indicated by a decrease in B(E2) values (see 
fig.[5|a). All three interactions suggest that the K 7T =5 + band in 90 Rh backbends at spin 
J m — 6=19 + . Backbending will be discussed after comparison of the calculated occupation 
numbers and E2 properties. 

The calculated j-shell occupancies (average number divided by the maximum number 
of particles in the shell) in 90 Rh are shown in fig||. To get particle occupancies in the SSSV 
calculation, the average number of holes was converted to the average number of particles. 
As can be seen from the figure, the p\/ 2 shell is filled almost all the time (occupancy is 
greater than 90% for most spin values) in the GF as well as SLGTO calculations, resulting 
in a half-filled gg/ 2 shell for all band states (except for the spin I m state in the GF 
calculation, as mentioned above). While in the SSSV calculation, this shell contains more 
than 10 particles (is more than half-filled) for spins J<21. The mixing reduces for 7>21, 
when 6 particles aligned in the gg/ 2 shell may contribute 21H. Note that an increase of 
the /V/2 shell occupation by particles, approaching the aligned state I m —16 in 48 Cr flUfl, 
corresponds to an increase of the gg/ 2 shell occupation by holes as spin increases to J m =25 
in 90 Rh in the SSSV calculation. 




Figure 3: Energies of odd-/ members of K™=5 + , T to t=0 band in 90 Rh calculated with 
GF, SLGTO (cf. [|l9j ) and SSSV interactions. For explanations see fig.||. Yrast states are 
connected by lines, while second excited states for each spin value are indicated by a 
corresponding symbol not connected by a line, except for the SSSV results shown in plot 
b, where the states forming bands (having large E2 transition strengths) are connected 
by dotted lines. 



The calculated E2 transition strengths in 90 Rh are shown in fig.[5]a. The B(E2) values 
calculated using SSSV interaction are larger than those obtained in GF or SLGTO calcu- 
lations, implying larger deformation and thus greater collectivity, as also seen from the 
energy plot with a rotational reference subtracted, fig.[3^. Additional reason for smaller E2 
transition strengths in more restricted calculations is that pi/ 2 does not contribute, thus 
there is only contribution from the g$/ 2 shell. In the case of half-filled j-shells, there is an 
additional suppression of the E2 transition strengths due to so-called "the center-of-the- 



shell selection rule" pq| : in a single- j configuration the spectroscopic moment vanishes, 
and the E2 transition can only take place between states differing in seniority by two 
units, if this j-shell is half-filled. Indeed, the occupancies (fig.|) show that the g 9 / 2 shell 
is almost always half-filled, Q spec ~ 0, and B(E2) are relatively small (figf|) in GF as well 
as in SLGTO calculation. The negative sign of Q sp ec calculated with the SSSV interaction 
(fig.||5) suggests prolate deformation along the band. 

All three calculations of 90 Rh predict a backbending at spin J m — 6=19 + , as would be 
expected from the discussion of band properties in a single- j shell assuming only JT=01 
pairing interaction being important (subsection |2.1|) . In the SSSV case, the backbending 
may be explained in terms of a band-crossing. The bands are shown in fig.^ where 
the calculated states are connected if there is a strong stretched E2 transition between 
them. As can be seen from the figure, yrare band crosses yrast band ( U K 7T =5 + band") 
between spins 7=19 and 21. For the spins below 21h, the excited band has lower particle 
occupation of the gg/ 2 shell (approximately 0.55) as compared to the yrast band value, 
~0.68 (fig.§). States 2lf 2 have almost identical occupation numbers. For higher spins 
the (79/2 shell occupation drops down to 0.5 in the yrast band, but equals 0.67 and 0.60 
at states with spins 232 and 252, respectively. The B(E2) values along the band with 
a lower occupation of the g$/ 2 shell are very similar to the GF and SLGTO results (cf 
fig.[|a), except for the transition 13^^11^ that is stronger, ~14.4 W.u. The Q sp ec values 
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Figure 4: Calculated shell occupancies in K=5 + band in 90 Rh. Three residual interactions 
were used: GF (marked by stars) as well as SLGTO (open circles) in the configuration 
space (pi/299/2) , and SSSV (filled circles) in (<79/2Pi/2P3/2) -10 - The j-shell is indicated 
close to a corresponding line. In the SSSV case, the number of holes was converted to the 
number of particles in shells. 



are approximately constant along this band, —10 e-fm 2 , while at 1=21 Q spe c decreases 
to —38.0 e-fm 2 . The continuation of "if=5 + band" has larger intra-band B(E2) values, 
than those calculated for the yrast states, except for the transition starting at spin 25^. 
The values are 1.4, 10.6 and 12.9 W.u. for transitions from state 25^ to state 23^, from 
23^ to 21^, and from 21^ to 19^, respectively. The spectroscopic quadrupole moment for 
the states J=21 2 -252 is varying around the value of —50 e-fm 2 . 

The B(E2) value calculated using SSSV interaction is reduced for the transition 
7i— 2), where I i =I m —A=21h. A straightforward explanation could be the fact that 
this transition connects two "rotational bands" having somewhat different structure, seen 
either in terms of the g 9 / 2 shell occupancy or the spectroscopic quadrupole moment. An 
alternative explanation could be the "center-of-the-shell" selection rule, mentioned above: 
the vanishing E2 transition would indicate purity of the half-filled single- j configuration 
[26|| . In subsection 2T it was discussed that the seniorities of states I m — 6 and I m — 4 are 
equal, thus the conditions for the selection rule are favourable. However, the Q spe c does 
not vanish neither at 19+ nor at 2lf (see fig.^|&), suggesting that the selection rule is not 
involved. While the structure of nuclear wavefunctions of those states is very different, 
and that is the reason for the reduced B(E2) value. 

A comparison of the predicted spin dependencies of the B(E2) values in 90 Rh (fig||) 
yields another question: whether the reduction of B(E2) at Ii=2lf is a feature of the 
SSSV interaction. The more restricted calculations do not predict reduction. On the 
other hand, they do not predict bandcrossing as a reason for the backbending. A sim- 
ilar situation occurs in 48 Cr calculated using the modified surface-delta interaction [12 



(MSDI) in the configuration space (/7/2P3/2) 8 , see fig.^|a. This interaction, fitted for this 
nucleus, describes excitation energies of the yrast states in 48 Cr a bit better than the KB3 
interaction does. The lowest bands predicted in 48 Cr by both KB3 and MSDI interactions 
are somewhat similar to those predicted in 90 Rh by the SSSV: In 48 Cr the yrast band 
continues above the backbending preserving the nuclear structure; however, the band, 



s 




Figure 5: Calculated E2 properties along the K 7T =5 + band in 90 Rh and the ground state 
band in 48 Cr: strengths of stretched E2 transitions, B(E2), are shown on the left; and 
spectroscopic moments, Q sp ec , are shown on the right. Values for 90 Rh were calculated us- 
ing the GF (marked by stars), SLGTO (open circles) and SSSV (filled circles) interactions. 
The modified surface-delta interaction (referred to as MSDI and marked by triangles) in 
the model space (/7/2P3/2) HU was used for 48 Cr, and the spin values for 48 Cr are shifted 
by 9h. 



terminating at spin J m =16 + , cannot be traced back below the backbending. The B(E2) 
values calculated in 48 Cr using MSDI are systematically smaller, most probably due to 
a restricted model space. Nevertheless there is a reduction of the B(E2) value for the 
transition starting at spin I{=I m — 4=12 in the MSDI calculation. (It is not clear why 
the Qspec vanishes at spin 12fr, see fig.[5]&. The selection rule would suggest, that it is an 
indication of a rather pure (iV/2) 8 configuration, however, the calculated probability of 
this configuration at this spin value is only 69%. But on the contrary to the selection rule 
prediction, the increasing probability of (iV/2) 8 is followed by the decrease of Qspec for 
spins 14 and 16.) However, the relative decrease of B(E2; 12^ — > 10^) was not observed 



\\22\\ . This decrease is also not predicted by the KB3 p4|| . Thus, having confronted the 
different predictions for the 90 Rh and 48 Cr nuclei, a conclusion may be drawn that only 
lifetime measurements or the calculation in a more extended model space may suggest a 
solution to the puzzle whether there is a decrease in B(E2; I{ — > h — 2) values at Ji=21 
in 90 Rh. 



2.3 Estimation of the role of pairing interactions 

The role of normalized L=0 pairing interactions |F] is discussed in present paper. Nota- 
tions for the forces used in the latter reference are adopted. Thus P01 ( or simply 01 as 
an index, e.g. Eqi) will stand for the normalized isovector J=0, S=0, T=l pairing inter- 
action, and -P10 (or 10) will stand for the normalized isoscalar J=l, S=l, T=0 pairing 
interaction. The pairing strength parameters are taken from eq.(C3) in ref. f7j. They are 
defined in terms of a constant multiplied by a scaling parameter, Tioj : Gqi=— 0.32Huj and 
Gio=— 0.51huj. An expression for Tiuj is taken from p^j , and the resulting strengths are 
given in Table |l[ 
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Table 1: Strength parameters of nor- 
malized pairing interactions used in 
present paper 



Pairing type 


22 Na 


48 Cr 


90 Rh 


Pqi-i Gqi 


-3.80 


-3.19 


-2.68 




-6.05 


-5.09 


-4.27 



The pairing interaction influence (i) to the state energy, (ii) the nuclear wavefunction 
and (iii) quadrupole properties will be discussed. Investigation is performed in a perturba- 
tive way, by comparing results obtained with the full Hamiltonian and the Hamiltonians 
having reduced content of some pairing interaction (either Poi or Pio). It will be shown 
that in some cases, a competition between the pairing and the rest of the residual inter- 
action may be seen. (Cf. discussion about the interplay between the quadrupole and the 
monopole parts of the Hamiltonian in ref . [f28[| . ) 

The Hamiltonian without Pjt pairing, Hjt(0), is defined as a full Hamiltonian, H, 
from which the normalized pairing interaction, Wjt, has been subtracted (in a similar 
way KB3-P01 and KB3-P10 were defined in ref. M, but the pairing interactions were not 
normalized. Note, however, that normalized and unnormalized forms are equivalent for a 
single major shell): 

H JT (0) = H-Wj T , (1) 

where JT=01 for the isovector (Poi) an d 10 for the isoscalar (Pio) L=0 pairing interaction; 
and the meaning of in parentheses will be seen later. To investigate the role of a JT 
pairing interaction in a perturbative way, lets define a Hamiltonian Hjt{X) that has some 
content of the pairing interaction indicated by a weight parameter A (0<A<1): 

Hj T (X) = H JT (0) + XWj T . (2) 

Thus Hjx{\) (=Hjt) is a linear interpolation between the Hamiltonian without pairing, 
Hjt{0), and the full Hamiltonian, Hjt{1)=H . One can parameterize the eigenvalues of 
this Hamiltonian in the following way: 

Ej T (X) = E JT (0) + (W JT )X + X (A), (3) 

where Ejt(0) is an eigenvalue of the Hamiltonian without Pjt pairing interaction, Ejt(X) 
is an eigenvalue of the Hamiltonian Hjt(X), (Wjt) is the expectation value of the normal- 
ized JT pairing interaction, and xW describes additional terms which are non-linear in 
A. One can see that x(0)=0. If x(A) is set to 0, the equation (^) describes first-order per- 
turbation to the energy Ejt{0). The energies are calculated and parameterized for each 
spin value separately, i.e. the quantities defined by eq.(H) have additional index - spin / 
- which is omitted (thus the eigenvalues in full notation are Ejt{X; I)). The expectation 
value of the pairing interaction, (Wjt), may be evaluated directly from the calculated 
nuclear wavefunction (cf. ref. Q). However, it was not possible to perform such evaluation 
in present investigation. 

An example of the Hjt eigenvalue spectrum dependence on the pairing interaction 
content is given in fig.|5|, where energies of the yrast states in 48 Cr calculated using the 
KB3 interaction are shown as a function of parameter A. The case when Poi interaction 
content is varied is shown on the left, while the Pio case is shown on the right. The 
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Figure 6: Calculated energies of yrast states in 48 Cr as a function of the Pjt pairing 
interaction content. KB3 interaction is used. The spin is indicated by a number close to 
the corresponding line. Plot a {b) shows the case JT=01 (10). 

meaning of parameters in eq.([3[) may be seen from this figure: (Wjt) is the slope of 
the Ejt{\) line at A=0, and xW is t ne deviation of Ejt(X) values from the linear 
behaviour (Ej T (0) + (Wjt)X)- One can see that |x(A)| increases together with A (the 
lines bend down), and that non-linear contributions are largest for lower spins. Since the 
slope of the E JT {\) line is proportional to the average number of JT-coupled pairs, it is 
clearly seen that there are more such pairs in the nuclear wavefunction calculated with 
a Hamiltonian close to the full Hamiltonian. Which may be understood from the fact 
that pairing interactions increase the binding energy of a pair, thus pairing correlations 
are more favourable if pairing force is stronger, i.e. the pairing interaction content in the 
Hamiltonian is larger. 

Lets define the pairing energy, AEjt, which will show the increase in the binding 
energy of a state having spin / due to an addition of the JT pairing interaction: 

AE JT = Ej T (0)-E(l), (4) 

where E(l) is the state energy calculated with the full Hamiltonian, and Ejt(0) is the 
state energy calculated using the Hamiltonian without Pjt pairing interaction; additional 
index - spin / - is omitted. The pairing energy will always be non-negative, since the 
pairing interactions are attractive, i.e. they increase binding. 

A Hamiltonian may be considered as consisting of monopole and multipole parts M, 
where the monopole part is responsible for the correct nuclear binding energies at the 
observed nuclear radii and for the single-particle behaviour, and the multipole (or resid- 
ual) part describes collective properties and is responsible for configuration mixing. One 
may notice that the pairing interactions contribute to both parts: the changed content of 
the pairing interaction affects the monopole part due to the moved centroids of the total 
interaction, while the multipole part is affected because this interaction is responsible 
for correlations (only this part of a pairing interaction, that contributes to the multipole 
Hamiltonian, may be viewed as being responsible for pair correlations in the spherical 
mean-field sense). The strengths of pairing interactions were derived in ref.0 from a 
monopole-free Hamiltonian. However, in present investigation it was not possible to make 
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exact separation of the pairing effects into monopole and multipole parts in a consistent 
way, thus only total effect is investigated. One term in the monopole part of the Hamil- 
tonian, namely, the single-particle contribution, may be readily evaluated having j-shell 
occupation numbers. The single-particle term may be written as: 

K d = Y.n 3 e v (5) 
{i} 

where Ej is the single-particle (or single-hole) energy of a j-shell, and rij is the average 
number of particles in the shell. 

In a similar manner as the change of binding energies (the pairing energy) is defined 
in eq.(f|), one can define the change of the single-particle contribution as an outcome of 
pairing interaction: 

AK d JT = K d JT (0)-K d (l), (6) 

where Kj T (0) and K d (l) is the single-particle term evaluated using occupation numbers 
calculated with the Hamiltonian without Pjt pairing and with the full Hamiltonian, re- 
spectively. This difference may give some idea about the change in the mixing with the 
high-lying j-shells as a result of pairing interaction, because particles occupying shells 
with large single-particle energy Sj make larger influence to the total single-particle con- 
tribution, see eq.(||). (For example, change in the occupation of the / 5 / 2 shell by 0.1 
particle and corresponding change in the fj/2 shell occupation yields 0.6 MeV in AK d , 
if the KB3 interaction is used.) In addition, it gives an idea of a fraction of the pairing 
energy coming due to single-particle contributions. 

The difference AK d T can be both positive and negative. If it is negative, more par- 
ticles are excited to higher j-shells in the full calculation, while the pairing interaction 
(over) compensates the loss in the single particle contribution. If it is positive, the pairing 
interaction contributes to the gain of energy by favouring configuration with less particles 
excited from the lowest j-shells. 



3 Results and discussion 



It is known that the eigenfunctions of a well-bound states are less sensitive to a slight 
change of the potential, than the corresponding eigenvalues are. If two shell model Hamil- 
tonians have two-body interaction being different by content of a pairing force, their 
eigenvalues are different but nuclear properties calculated using the corresponding eigen- 
functions may be similar. The discussion of the role of pairing interaction is started from 
investigation of this interaction's influence on the nuclear wavefunction and quadrupole 
properties (subsection |3.1|). It is known, that the pairing interaction does not affect 



quadrupole properties [2S, 29 1. It is, however, shown here that there are cases when this 
is not true. The influence of pairing interactions on energies is discussed afterwards, in 



subsection |3]2|. A strict separation into these two topics cannot be made, thus the energies 
are also discussed together with wavef unctions. 
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Figure 7: Energies, occupancies of the shells, single-particle contributions, and E2 prop- 
erties of yrast and yrare states in 48 Cr calculated at spins J 7r =0 + , 4 + and 6 + using the 
KB3 interaction. The Poi pairing interaction content is given on the horizontal axis 
(A=0 means no-pairing interaction, while A=l means full pairing interaction is present 
in the Hamiltonian). The top- row plots show the energies, E; the second-row plots show 
shell occupancies; the third-row plots show the single-particle contributions, Kq ± ] and the 
bottom-row plots show E2 properties: E2 transition strength between yrast 2 + and + 
states, and the spectroscopic quadrupole moment Q S pcc at spins J=4 + and 6 + . The solid 
and dot-dashed lines are used to distinguish the yrast and yrare states, respectively. The 
f 7 / 2 shell occupancy is marked by diamonds, p 3 / 2 by crosses, pi/ 2 by pluses, and f 5 / 2 by 
triangles. 
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3.1 Nuclear wavefunction and quadrupole properties 

In general, the change in energy of a state Ejt(X; I) as pairing content (i.e. the A value) 
increases, illustrates two effects. Besides a direct effect coming from the increased binding 
of certain nucleon pairs (linear in A), there is an additional increase in energy (non-linear) 
due to the change of the number of such pairs that also results in abandoning of some less 
favourable configurations, i.e. the "structural effect" or x m ec L-(|!)- From fig.|6|one may see 
that the structural effect increases the binding energy: if we start from the Hamiltonian 
without, say, Poi pairing, there are not so many JT =01 coupled pairs, and the slope is not 
large. But as A approaches 1, it becomes more and more favourable to form such pairs, 
and the P(A; I) curve bends down. If we start from the full Hamiltonian and reduce the 
amount of Poi pairing, the energy loss would be much larger if the nuclear wave function 
is fixed than it is when configuration change is allowed. 

Investigation of the pairing interaction effect on a particular state in a perturbative 
way may reveal a gradual change of the state properties. Changes are larger, when there is 
another close lying state, since these states may mix in a different proportion following the 
change of the interaction. As a result, when varying A, one can see a mutual interchange 
of state properties. In an extreme situation, either sharp or gradual transition of the state 
properties from those of one state to those of the other state may be observed. 

Let us first investigate a case, when there is a substantial energy distance between two 
states, thus other states do not influence the change of the state properties. The fig.|7| 
shows some properties of the first and second excited states in 48 Cr calculated at spins 1=0, 
4 and 6 + using the KB3 interaction, and the content of Poi interaction is varied. Those 
properties are: state energy, j-shell occupation numbers, single-particle contribution to 
the total energy, and spectroscopic quadrupole moment Q spe c (or E2 transition strength 
where Q spe c is not defined). It is seen that the energy of the yrast states decreases (binding 
increases) as A increases, but the occupancies almost do not change (except for the p 3 / 2 
shell occupancy at + ). Some occupancies are difficult to distinguish in the plot, however, 
their values are not important for the discussion. The occupancies may be summarized 
in terms of the single particle contribution to the state energy, see eq.(|). As can be seen 
from the figure, of the yrast states changes only slightly or not at all. One may also 
notice that the wavefunction of the yrare state at + is more affected by the change in 
pairing force than the yrast state (occupation numbers and thus Kq X change more than 
in the yrast state). In addition, the second state at 7=4 + interacts with the third state, 
and the occupancies, the single-particle contribution as well as spectroscopic quadrupole 
moment change drastically at A~0.45, while the state energy only slightly changes slope. 
A similar change is seen at 6^. 

The bottom-row plots in fig.|7| show calculated E2 properties in 48 Cr as a function of Poi 
pairing interaction content. As may be seen, the behaviour of the spectroscopic moment, 
Qspec, is very similar to that of the single-particle contribution, Kq X : if does not 
change, the Q sp ec does the same. However, it is wrong to correlate the increasing f 7 / 2 shell 
occupancy (K^ tends to 0) with the more oblate nuclear shape (Q sp ec>0). The "center-of- 
the-shell" rule (mentioned in subsection above) demands Q spe c=0 for K d =0, however, 
as can be seen from the figure, small admixtures of other j-shells causes deviation from 
the rule. The bottom-left plot shows the strength of the stretched E2 transition between 
the lowest yrast states, 2f— >0f. As more pairing interaction is present in the Hamiltonian 
(A increases), less particles are excited from the fj/2 shell (K^ decreases), thus prolate 
nuclear deformation decreases (cf. the behaviour of Kq X an d Qspec at 1=4%) ■ This is in 
agreement with the "center-of-the-shell" rule. However, the low probability of pure f~, 2 
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Figure 8: Calculated probabilities of configurations in the 48 Cr ground state as a function 
of Pjt pairing interaction content, A: the left and right plots show P i an d Plo cases, 
respectively. The KB3 interaction was used in the calculation. Only configurations with 
probability exceeding 5% at some point are shown. The missing probability is about 33%. 

configuration (see below) as well as changing Q spC c of 2f state (not shown here) suggests 
that a bit more sophisticated picture than the mentioned selection rule is needed. Note, 
however, that the reduction of the B(E2) value by 6% means only insignificant change in 
an approximate quadrupole deformation j3 2 (calculated assuming the axially symmetric 
rotor model is valid): from ,$2=0. 270 to 0.263. 

Although the binding energy of the yrast state changes significantly due to a change 
in the pairing interaction content, the occupation numbers change only slightly, and only 
a more detailed view of the probabilities of different configurations reveal changes. As 
an example, configurations in the ground state of 48 Cr are shown in fig||. Those con- 
figurations which never achieved probability larger than 5% are omitted. The missing 
probability, i.e. summed probability of not shown configurations, constitutes 35% and 31- 
35% for P 01 and P w cases, respectively. The average numbers of particles in j-shells are 
unevenly distributed among the omitted configurations. For example, the missing particle 
numbers are 1.8 for f 7 / 2 , 0.5 for p 3 / 2 , 0.2 for pi/ 2 and 0.4 for f 5 / 2 shell in the case of P i 
pairing interaction. The probabilities marked like fy 2 ®fy 2 P 3 / 2 are added probabilities 
of configurations ^fy 2 ^{.fy 2 P 3 / 2 ) and ^{.fy 2 P-3,/ 2 )^fy 2 (here n stands for protons and v 
for neutrons), while marks like {fy 2 ) 2 indicate that protons and neutrons have the same 
configuration {fy 2 in this example). The fig.||a shows that the probability of configura- 
tions where all 8 particles are in the fi/ 2 shell is increasing by 8% in the Pqi case, when 
A goes from (no pairing) to 1 (full pairing). At the same time, the configurations with 
one particle excited to p 3 / 2 decreases by the same amount. The P i pairing increases the 
probability of configurations with two particles of the same kind being excited to f$/ 2 by 
3%, but this probability reaches only about 6% at the maximum. The mixing with p%/ 2 
may be favourable at A=0, but at A=l, where JT=01 pairs are more favoured in energy, 
shells with large £, i.e. f-j/ 2 and f$/ 2 , are more important. The fig||& shows that the Pio 
interaction lowers the probability of fy 2 as well as fy 2 p 3 / 2 configurations by 4% each; 
while the probabilities of configurations with one or two particles in the f^/ 2 shell slightly 
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increase. 

Information about changes in the nuclear wavefunction are still hidden in fig.[8[ the 
linear change of configuration probabilities does not explain non-linear change of the state 
energy (see fig.|6]). It may be suspected that the average seniority changes in a nonlinear 
manner. 

An example, where the two lowest-lying states affect each other when the pairing in- 
teraction content changes, may be taken from 90 Rh calculated using the SSSV interaction. 
Fig.p illustrates three cases of the two interacting states calculated for the Poi force: weak 
mutual influence at spin P r =13 + , smooth "crossing" at 17 + , and sharp crossing at 19 + . 
The states hardly influence each other at spin J=13 + (see left column plots in figH). Al- 
though the energy difference between the states remains approximately constant as more 
pairing is added (i.e., A approaches 1), properties of the states change in a correlated 
way: in the yrast state the #9/2 shell becomes more mixed with other j-shells, while the 
opposite is seen in the yrare state. This change in shell mixing can be seen both from 
the shell occupancies by particles and from the single-hole contribution, K^. Note that 
the single-hole contribution in 90 Rh behaves in an opposite way than the single-particle 
contribution in 48 Cr when trend is compared to the occupation of the high-j shell: if the 
fr/2 shell occupancy increases, cf. fig.0, the single-particle contribution in 48 Cr decreases, 
while the gg/2 occupancy by particles and the single-hole contribution in 90 Rh increase or 
decrease together. The two lowest states come closer in energy and separate again at spin 
I=17 + (see middle column plots in fig.^j). The behaviour of occupation numbers (or of 
the single-hole contribution) suggests that the states gradually exchange their properties, 
since the shell occupancies of one state at A=0 become equal to those of the other state 
at A=l. A similar picture is seen at spin 7=19 + where the two states approach and cross 
each other without much interaction: the energy slope changes sharply at A=0.9 and both 
the occupancies and thus the single-hole contributions change abruptly. In this case it is 
possible to calculate the pairing energy of a specific configuration, i.e. of the configuration 
being yrast at A=l, or of the one being yrast at A=0. 

The bottom-row plots of fig.|9| show Q spe c calculated in 90 Rh using the SSSV interac- 
tion for the three spin values discussed above. It is clearly seen that Q sp ec is correlated 
with the single-hole contribution, Kq V For the less mixed configurations (small 
values), Q S p e c tends to zero; while for the more mixed configurations (larger K^), the 
prolate deformation increases (Q spe c decreases). A picture of the mutual influence be- 
tween the two states, obtained plotting the deformation coordinate (or Q spec , assuming 
the axial symmetry) at horizontal axis and energy at vertical axis, would be somewhat 
different than the one described above. If the positions of points representing values of 
first and second states for a given spin are followed as A goes from to 1, one would 
see, that the "minima" are well separated by deformation, and they merely "sink" (i.e., 
the binding increases) at spin 13 + . A bit more complicated picture occurs at spin 17 + , 
where minima are separated for A values close to the limits, but at A~0.7 the two states 
should be treated as first and second state in the same potential well and not as sep- 
arate minima. The properties of states are also very similar there (see fig|| middle 
column). At spin 19 + , there minima are well separated by deformation for A=0. But 
for larger A values, the well depths increase at a different rate and become more and 
more similar. As a result, the wavefunctions start to mix, and finally, the state having 
larger deformation becomes yrast. A similar sharp transition between two states may 
occur in the backbending regions, if the crossing bands have different pairing properties. 
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Figure 9: Some properties of yrast and yrare states in Rh calculated at spins 7=13 , 
17 + and 19 + , using the SSSV interaction, as a function of the Poi pairing interaction 
content. For explanations see fig.|7|. The particle occupancy of the gg/2 shell is shown by 
filled circles, by crosses, py2 by pluses. The single-hole contribution is shown in 
the third-row plots. 

From the presented here cases of two interacting states one may see, that the ( JT=01) 
pairing increases shell mixing and may unfavour the pure high-? shell configurations. Thus 
deviations from the single-? picture of the pairing energy spin dependence, discussed in 



subsection [2J] , may be expected. Although only the Poi cases were discussed here, a 
similar picture is obtained when P 10 interaction content is varied. 



17 



3.2 State energy 



In the previous subsection it was already discussed that the increase of pairing interaction 
content yields additional binding of the state. This subsection discusses changes in the 
binding energy of yrast states. First of all, investigation of the obtained pairing energy 
dependence on the model space is presented. Then pairing energy dependence on the 
residual interaction is briefly reviewed. And finally, the pairing energy in the selected 
bands is discussed. 



3.2.1 Model space effects 

During discussion of the role of pairing interactions, a natural question may be asked, 
whether the model space is sufficient to trust the results. For example, the definition of 
schematic pairing forces {L=0) suggests that it is important to include both the high-j 
shell and its spin-orbit partner in the model space, in order to get reasonable pairing 
energy values. Thus the model space for 90 Rh is probably too small, since the #7/2 shell 
is excluded. This section tries to shed some light on this problem. 

To investigate the question whether the pairing energy, defined by eq.(Q), is model 
space independent, yrast states of 48 Cr were calculated in different model spaces made 
of different combinations of j-shells from the pf major shell but using the same residual 
interaction, KB3. Experimental energies of the yrast states are reproduced with a different 
success in the selected model spaces, but in this investigation only the response to the 
change of pairing interaction content is important. The derived pairing energies (defined 
as energy difference between yrast states calculated with the Hamiltonian having no Pjt 
pairing interaction and with the full Hamiltonian) are shown in fig.|T0|: Pqi case is shown 



on the left, and P\o case is shown on the right. It may be seen that the largest contribution 
comes from the spin-orbit partners having large orbital angular momentum I. Thus the 
change in binding energies calculated in the full pf model space is almost completely 
accounted for in a smaller model space containing both and f^/2- It is somewhat 
surprising that there are no obvious signs of a model space increase from /V/2,/5/2 to 
full pf for both P i and P10 forces. This may be understood taking into account the 
relatively small orbital angular momentum of the p%/2 and pi/2 shells, as well as from 
small occupancies of the j-shells other than the high-j shell (see fig.^and ref.||14||). 

The similarity of the pairing energy in the /V/2,/5/2 and full pf calculations suggests 
that there is some saturation of the pairing interaction influence, and the change in nucleus 
structure accounts for the interaction change. That is also suggested by the model space 
dependence of the change in the single particle contributions (see below). 

One conclusion which could be drawn from this investigation is that the approximate 
energy change due to the P01 force as a function of spin can be obtained from the calcula- 
tion in a model space without the spin-orbit partner, for example, /7/2P3/2 or /7/2P3/2P1/2 
(see fig.|T0|a) . The absolute magnitude can be estimated multiplying the pairing energy 
obtained in a more restricted model space by the ratio of model space degeneracies. For 
example, if one takes the result of the /7/2P3/2P1/2 model space and multiplies it by 20/14, 
the estimated AEqi for the full pf calculation is only slightly too large: the difference 
is about 0.3 MeV for spins 1=0 and 2 + and less for other spin values. However, the 
predicted slope between states having spins 1=8 and 10 is different from that in the full 
calculation. Such scaling does not work for the P\q force (see fig. [[(]&), where inclusion 
of f§/2 is necessary to obtain the profile of the AE W spin-dependence. One may also 
notice, that the inclusion of the f 5 / 2 shell into the model space approximately doubled 
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Figure 10: Model space dependence of the pairing energy. Calculated for the 48 Cr nucleus 
using KB3 interaction. The left (right) plot shows Poi (-Po) case. 

the calculated P w pairing energy. 

This investigation, however, leaves the question, whether inclusion of the other major 
shells would change the pairing energy. One would expect that the influence of higher lying 
j-shells is not large because a large difference in single-particle energies would limit the 
mixing. If, however, the shells below the lowest partially-filled shell (the "lower shells") 
are included in the model space, the pairing interaction effect should be separated into a 
part coming from the lower shells and the rest, because the lower shells always react to 
the change of the pairing interaction content (since they are mostly filled). 

The model space dependence of a difference in the single-particle contribution, AKj T , 
defined by eq.(^|), is shown in fig .[11]. This difference gives some idea about the change 
in the mixing with the high-lying j-shells as a result of pairing interaction. It can be 
seen that AKj T behaves in a rather irregular way, as compared to the total change in 
energy (fig.[L0]), and it is not possible to predict how the contribution will change with 
an increase of the model space. However, there are also similarities in the change of the 
single-particle contributions. For example, it is clearly seen, that the pairing interactions 
favour configurations with more particles excited to higher shells in the model spaces 
containing both / 7 / 2 and / 5 / 2 shells (except for I=0 + and 12 + states calculated for the 
Poi case in the full pf model space), since AKj T <0. The opposite is seen for most spin 
values in the /7/2P3/2 -Poi case. It is also interesting to note that the P w force decreases 
shell mixing at the maximum spin state 7 m =16 + , since AKj T is positive here. 

In summary, the investigation of the model space dependence of the pairing energy 
suggests that the model space containing shells slightly above the lowest partially-filled 
shell is sufficient to investigate the influence of the Poi pairing interaction on the state 
energies, while the P\q pairing energy can be obtained only in model spaces containing 
both the high-j shell and its spin-orbit partner. 



3.2.2 Residual interaction dependence 

Somewhat different investigation is the calculation of energy changes in the same nucleus 
using different interactions. If the residual force is unique, the difference in the calculated 
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Figure 11: Model space dependence of the difference in single-particle contribution, 
AKj T =Kj T (0)-Kj T (l). Calculated for nucleus 48 Cr using the KB3 interaction. The 
left (right) plot shows Pqi (Ao) case. 



pairing energy would come only from the model space difference. However, since there is 
no unique force yet, the agreement in the pairing energy, if occurs, is more likely to be 
accidental, and there should be clear dependence on the chosen residual interaction. A 
separate question is whether this dependence should show up only as different nonlinear 
effects, x m ec i-(3)- If t ne nuclear mean-field (Hartree-Fock field) is universal (cf. discussion 
about the universal multipole part of the Hamiltonian in ref.[[7|, that gave foundation for 
this investigation), one would expect that the average number of pairs in the wavefunction 
calculated using the Hamiltonian without pairing interaction is independent on the chosen 
(full) residual interaction. A disagreement in this average number of pairs would come 
from the other multipole terms favouring certain two-nucleon states, different in different 
residual interactions. The average number of pairs is proportional to (Wjt) in eq.(|3]), 
thus more information will be provided in the subsection |3.2.3| , when the coefficients of 
the state energy decomposition, eq.(|3]), will be presented. 

The investigation was performed for 90 Rh using GF and SLGTO interactions in the 
configuration space (pi/25'9/2) 14 , as well as the SSSV interaction in (g9/2Pi/2P3/2)~ 10 - The 
calculated pairing energies are shown in fig.[12|. To take into account the effect of the lower 
shells (in our case, of thepi/2 shell), which was discussed in subsection |3.2.1| , pairing energy 
calculated in 90 Rh using GF and SLGTO interactions was approximately corrected taking 
the maximum spin state 7 m =25 + as a reference. There spin is produced by 10 particles 
in g 9 /2 shell, and the pi/% shell is full. There is no non-linear dependence on the pairing 
interaction content, i.e. x(A)=0 at this state. Thus the values obtained in the GF and 
SLGTO calculations are shifted down to match the SSSV value at the maximum spin 
state J m =25 + . GF pairing energies are reduced by 0.66 and 0.24 MeV for P i an d -Pio 
cases respectively. The corresponding values for the SLGTO interaction are 0.80 and 0.12 
MeV. As explained in the subsection [3.1| above, due to a sharp "crossing" in the energies 
Ejt(X',I) calculated with the SSSV interaction, it is possible to calculate the pairing 
energy of configurations becoming yrast in the full Hamiltonian and being yrare in the 
Hamiltonian without pairing. These values, ^(O)— Ei(l) (here subscript refers to the 
state number for a given spin and total isospin), calculated at spins J=15 + and 19 + are 
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Figure 12: Interaction dependence of the pairing energy calculated for 90 Rh nucleus using 
three interactions. The left (right) plot shows Poi (Plo) case. The GF and SLGTO values 
are shifted down to match the SSSV value at spin J=25 + . For spins J=15 + and 19 + , the 
SSSV energy differences between yrare at A=0 and yrast at A=l states are indicated by 
cross signs, and the energy differences between states behaving in an opposite way are 
marked by plus signs. 

marked by crosses in the fig.0. The pairing energy of configurations which become yrare 
in the full Hamiltonian, Ei(0)—E 2 (l), are shown by pluses. The figure clearly shows that 
the configuration becoming yrast in the calculation with the full interaction gets much 
more pairing energy than the yrare configuration. 

As seen from the figure |12|a, the spin dependence of Poi pairing energies, obtained from 
different interactions, is very similar, especially between the GF and SLGTO interactions 
calculated using the same model space. However, the absolute magnitude of pairing 
energy is different. The difference in the calculated Poi pairing energies using GF and 
SSSV interactions cannot be explained by the model space effects (cf. fig JlOD : AP 01 from 
the SSSV calculation should be larger than values obtained from the GF calculation. Thus 
it is the expected dependence on the residual interaction. 

The origin of the difference in the spin dependence of Poi pairing energy close to the 
band termination is not clear: in the SSSV case, there is a small increase in pairing energy 
at spin 21 + (in a similar way as in 48 Cr at I m — 4=12, see fig. |lO|) , while in more restricted 
calculations using GF and SLGTO interactions, the Poi pairing energy is locally smaller 
at 21 + and increases at i=23 + . 



The above discussion of the model space effects (subsection |3.2.1| ) suggested, that the 
calculated Pio pairing energies should be smooth if model space does not contain the 
spin-orbit partner of the high-j shell. And indeed, the APio values calculated for 90 Rh 



with GF and SLGTO interactions are smooth (see fig.|i~2"|&) , like in model spaces without 
fs/2 shell for 48 Cr, see fig.|H]. The SSSV calculated energies behave in a different way, 
suggesting the anticipated interaction dependence of the calculated pairing energy 

There are more cases when different interactions predict similar pairing energy. For 
example, the Poi pairing energy calculated along the ground state band in 48 Cr using the 
modified surface-delta interaction in the model space (/7/2P1/2) UK is very similar to the 
values calculated using the KB3 interaction in the full pf model space. In addition, even 
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Figure 13: Interaction dependence of the single-particle contribution difference calculated 
for 90 Rh nucleus using GF, SLGTO and SSSV interactions. In the latter case the single- 
hole contribution difference is shown. The left (right) plot shows Pqi (Plo) case. The 



notation is the same as in fig. 12 



the decomposition into direct and structural effects (eq.@) gives almost identical values. 
The MSDI parameters were fitted to reproduce experimental energies in this nucleus, thus 
the model space effect, i.e. the absence of f 5 / 2 , was probably remedied by this adjustment. 
The Pio pairing energy calculated with MSDI is smaller than the one obtained using KB3 
in the same model space, shown in fig.|T0|. 

Fitting to the energies, in general, does not imply that all energy-related properties 
will be well described. For example, the Pqi pairing energy, calculated for 48 Cr using the 
surface-delta interaction fitted to the energies in the full pf model space (the parameter 
values At=0.2+0AT), is larger than KB3 values. However, the spin dependence has the 
same profile. 

The single-particle contribution difference, calculated in 90 Rh using three different 



residual interactions, is shown in fig.|T3|. As was already seen from the model space in- 
vestigation of AK d , the contributions may be different in these calculations, first of all, 
because of the different model spaces. Thus the K d values obtained in GF and SLGTO 
interactions, being one similar to another, should be expected to be different from the the 
SSSV values. As seen from the figure, the pairing interactions favour configurations with 
more holes excited from the gg/2 shell in the SSSV calculation. A similar conclusion was 
made for 48 Cr. 



3.2.3 Pairing energy 

The JT pairing energies, calculated for three selected nuclei, are shown in fig.fni The spin 
values in 22 Na and 48 Cr calculations are shifted to match the maximum spin value 7 m =25 + 
of the selected band in 90 Rh by 16 and 9h, respectively. Pairing energies in 22 Na and 48 Cr 
are shown as obtained from the calculations using USD and KB3 interactions, respectively. 
While for the nucleus 90 Rh, the estimated Poi pairing energy is shown. Estimation tries 
to account for the spin-orbit partner of the gg/2 shell, namely the g-?/2 shell, that is not 
included in the model space. It is made by assuming that the pairing energy scales with 
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Figure 14: Spin dependence of the calculated pairing energy in the selected bands of 22 Na, 
48 Cr and 90 Rh nuclei. For explanations see text. The Poi (Plo) case is shown on the left 
(right). 



the model space degeneracy, as suggested by the model space investigation in subsection 



3.2.1. The estimation is based on the calculated SSSV pairing energies (since they are 



obtained in a larger model space than, say, GF values), thus the scaling parameter is 1.5. 
The Pio pairing energy in 90 Rh is taken as obtained using the SSSV interaction, because 
Pio energies do not scale with the model space. As can be seen from the figure, both 
the P i an d Pio pairing energies decrease as spin increases (on the average). In addition, 
the values are similar in magnitude for the odd-odd nuclei, 22 Na and 90 Rh. In the 90 Rh 
calculation using the SSSV interaction AP 10 is even larger than AP i- While in the 
case of even-even nucleus, 48 Cr, the influence of JT=01 interaction is larger, as was also 
shown in ref. pf. Thus both pairing interactions reduce the moment of inertia. However, 
the relative importance of Poi and Pio pairings to the moment of inertia in N=Z nuclei 
depends on whether it is an odd-odd or even-even nucleus; in addition, it depends on the 
total isospin of states as well as on the oddity of state spin. Fig.|TJ] shows that the spin 
dependencies of Poi and Pio pairing energies are more or less correlated. Those issues are 
addressed in a greater detail in ref.[|ll|] where pairing properties of the rotational bands 
in the odd-odd N=Z nucleus 62 Ga are discussed. 

Due to a spin shift when plotting, it is easier to see that the irregularities in the spin 
dependence of the Poi pairing energies are similar in the bands (see fig.[T4]a). There is 
no JT=01 coupled pairs in the aligned, band terminating state, thus there is no pairing 
energy at spin J m . There is a local maximum at spin I m — 4, which happens to be the 
backbending point in 48 Cr and 90 Rh (see figs. and [3]). This state was also described 
as the maximum spin state for seniority f max — 4 (where t> ma x=2j+l), and the discussion 



in subsection |2.1| suggests that the pairing energy should be smaller in the neighbouring 
states: at higher spins the total seniority becomes larger; while at spin J m — 6 seniority is 
the same as in state J m — 4 in a single- j picture, however, this state is less resistant to the 
shell mixing that breaks JT=01 pairs, since it is not "maximum spin for a given seniority" 
state [0. The Poi pairing energy increases as spin decreases from J m — 6 to 7 m — 10, and 
continues to increase for 48 Cr towards the lower spins, but remains rather constant for 
90 Rh. It may be an effect from the T=l pairing blocking in T tot =0 states of an odd-odd 
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nucleus. 

The spin shift reveals also similarities of the AE 10 spin-behaviour when values are 
compared starting from the band terminating state, see fig.|i~4|&. The J m — 4 state has an 
irregularity as compared with the trend of lower spin states. For the spins below J m — 4, 
the AEiq is almost a linear function of spin. In 48 Cr there is a slight increase of AE\o (as 
compared with this linear interpolation) at spin I m — 12=4 + . This spin value is the lowest 
state if only P w pairing interaction constitutes the complete interaction. Thus it may 
be a signature, that KB3 contains overestimated P\o force (cf. fig.[l]), because a similar 
irregularity is also seen in the full energy spectrum with a rotational reference subtracted, 
while there is no such irregularity in experimental data (see fig-D). 

The parameter values in the decomposition of Pjt pairing energies, see eqs.(^) and 
(fD, are shown in fig .[15]: plots in the top row show expectation values of the normalized 
pairing interaction (Wjt), while the bottom row plots show the total structural effect 
Xjt(1)- Like in fig.[T4], the spin values for 22 Na and 48 Cr are shifted. The parameter values 
for 22 Na and 48 Cr are calculated using USD and KB3 interactions, respectively. For the 
90 Rh nucleus, parameter values obtained using all three interactions, GF, SLGTO as well 
as SSSV, are given. The values are not scaled to estimate the effect from the excluded #7/2 
shell. The irregularities of x(l) values calculated for 90 Rh may be attributed to the two 



state interaction, discussed in section |3.1| . To eliminate the effect coming from particles 
in the pi/2 shell, the (Wjt) values obtained in the GF as well as SLGTO calculations are 
shifted down to match the SSSV value at i=25 + (cf. fig.|i~2"D. The x(l) values show that 
the non-linear term, and thus the reconfiguration effect, is not negligible and plays an 
important role in shaping the spin dependence of the bands. In particular, there is an 
increase in the nonlinear contribution at spin (J m — 4)h. The difference in the magnitude of 
the parameter values in 22 Na, 48 Cr and 90 Rh cannot be explained by simple mass scaling, 
hw . There are additional effects, that may have origin in the dependence on the residual 
interaction. 

As seen from the fig.|T5]a, the expectation value of the P01 interaction is approximately 
constant for spins ranging from I m — 2 to I m — 6 (or even I m — 8), suggesting that the average 
number of JT=01 pairs in this spin range is not changing in the Hamiltonian without 
Poi pairing interaction. The fig.^c suggests, that the structural effect is responsible 
for the increase of the pairing energy at spin I m — 4. The figs.|T5]a and c suggest, that the 
unfavoured band termination is due to the equal number of JT=01 pairs at spins I=I m —A 
and J m — 2 in the calculation without this pairing interaction, since the structural effect 
changes by the same amount when going from J m — 4 to I m — 2, and finally I m state. The 
T=0 pairing helps in making the I m state unfavourable: here the change in the pairing 
energy is larger when going from I m — 2 to I m state as compared to the change between 
states I m — 4 and I m — 2. 

The structural change is also important for the effect coming from the Pio interaction 
(see plots in the right-column of fig. 13). The expectation value of the pairing interaction 
has smoother spin dependence than the pairing energy does in the odd-odd nuclei (cf. 
fig.|14|). It is interesting to note, that the difference in AEiq between the GF (or SLGTO) 
and SSSV interactions is coming from the structural effects, while the average number of 
L=0, S=J=1, T=0 pairs is the same at A=0. 
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4 Summary and conclusions 

There are several ambiguities, when discussing the role of the pairing interaction in the 
shell model. First problem is the "no-pairing" Hamiltonian. Since it is defined as a 
deviation from the full Hamiltonian (see eq.(|])), the correctness of defining the magnitude 
of this deviation (i.e. the pairing strength) is crucial to have a correct Hamiltonian without 
pairing. Having selected (or fixed from other considerations) the strength of the pairing 
interactions, a question may be asked whether a full force is present in a chosen residual 
interaction. This question was also discussed in ref.fTj, where it was argued that the 
multipole part of realistic interaction, that is responsible for particle correlations, including 
pairing interactions, is almost the same. Thus only small deviations from the pairing 
interaction content, estimated in this reference, are expected. However, the model space 
dependence of the pairing properties may yield unwanted effects. 

Another problem is the form of the pairing interaction. The normalized JT=01 pairing 
interaction was shown |7| to remedy the divergence in the pairing plus quadrupole model, 
in the case when many shells are included in the model space. However, this normalized 
form still does not solve the divergence problem completely. Thus the form of the force 
is not final yet, and details of the derived role of pairing interactions may change. 

The third problem is the interpretation of the role of a pairing interaction in compar- 
ison with that derived from other models. The difference in treatment of pairing corre- 
lations in the shell model and, say, Hartree-Fock-Bogolyubov approach puts a question 
whether comparison is meaningful. 

The fourth problem is the relation between the pairing properties obtained in the shell 



model and observables. 

Having the above ambiguities in mind, lets summarize the presented investigation of 
the role of pairing interactions. Three selected bands in 22 Na, 48 Cr and 90 Rh, having the 
total isospin T tot =0, were discussed. Their choice, based on the behaviour of the JT=01 
pairing energy in a single- j shell model, was explained. Since there is no experimental 
data on 90 Rh available, three interactions were used to calculate its properties, and a 
comparison of their results was presented. All three interactions predict that the ground 
state is 1=0, T to t=l, and that the excitation energy to the lowest T tot =0 states is ap- 
proximately 1.3 MeV. In addition, the K n =5 + (T tot =0) band, which was of interest in 
the present investigation, backbends at spin 7=19 + . In the SSSV calculation, the back- 
bending is related to a bandcrossing. The SSSV calculation should describe quadrupole 
properties more reliable, thus odd-/ members of the K=5 + band are expected to have 
a prolate deformation. The comparison of the energies calculated in 90 Rh using the GF 
(or SLGTO) and SSSV interactions raises a question, whether the model space used to 
describe 48 Cr yrast band is large enough to discuss also the high spin states. This ques- 
tion is supported by the observation that the deviation between calculated and observed 
energies is increasing after 1=10 towards the aligned state at 16 + (see fig||). In addition, 
the proposed shell model explanation for the backbending as a result of interplay between 
the monopole and quadrupole forces p8| does not exclude the picture of the unpaired 
band-crossing (ground state band crosses with the K n =8 + band). 

A perturbative approach to investigate the role of schematic normalized L=0 pairing 
interactions J7| was presented. Both T=0 and T=l pairing interactions reduce the moment 
of inertia (that is in agreement with ||), however their relative importance have interesting 
dependence on the spin oddity and total isospin of the states (see [0). It was shown 
that the pairing force brings energy in two ways: (i) directly, due to an extra binding of 
certain pairs of nucleons, and (ii) indirectly, due to an increased likelihood to form those 
pairs. From the calculations it follows, that the indirect effect is as important as the 
direct effect. 

Problems occurring during investigation of the role of pairing interactions were dis- 
cussed. In particular, the case of two close-lying states was described. It was shown that 
the pairing interaction do not alter quadrupole properties, if the energy distance between 
the states is large enough. 

The dependence of the pairing energy (the difference between energies of yrast states 
calculated from the Hamiltonian without pairing interaction and the full Hamiltonian) 
on the model space was investigated taking 48 Cr nucleus as an example. It was shown 
that the spin dependence of the Poi pairing energy derived in the full model space may 
be estimated from the calculation in a smaller model space by taking into account the 
increased degeneracy of the model space. This observation was used to estimate the 
JT=01 pairing energy in 90 Rh as would be obtained if the (77/2 shell would be included 
in the model space. It was also shown that the P10 pairing energy does not scale with 
the model space degeneracy. The pairing energy dependence on the choice of the residual 
interaction was also discussed. During the investigations of those two dependencies, it 
was shown that the pairing interactions in most cases facilitate the shell mixing (a similar 
picture occurs in the BCS model). 

The pairing energies calculated in the three selected bands were compared. It was 
shown, that if the spin values are shifted so that the aligned state having spin I m =(p+1) 2 
coincides, the spin dependencies of pairing energies in the three selected bands have similar 
profiles. The decomposition of the pairing energy into direct and indirect parts, available 
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due to a perturbative way of investigation, was also presented. It was shown that the 
similarities in the spin dependence of the coefficients remain. 

Despite being rather detailed, this investigation raised questions which could not be 
answered. For example, the energy spectrum of 90 Rh remains an open question. The 
three used interactions predict different not only excitation energies, but also the energy 
distance between the yrast and yrare states. The origin of the backbend of the K n =5 + 
band in this nucleus is also unclear: whether it is similar to the one discussed in ref. p8 



or it is "usual" backbending originating from a band-crossing. The investigation of the 
pairing energy dependence on the model space suggests that the calculation performed 
in the model space containing also #7/2 shell may give somewhat different results than 
discussed here. The experimental data would provide a stringent test for the choice of 
the residual interaction. 

In addition, a systematic study of the pairing energy dependence on the total isospin 
and the oddity of state spin is needed to add an understanding what is the origin of 
blocking effects. 
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